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1. Introduction 

It is well known that a comprehensive theory of strong interactions is still in an embryonic stage, as we 
understand little about non-perturbative QCD. On the other hand, considerable progress has been made 
in building models for soft scattering at high energies [1-5]. The main constituent of these models is 
the soft Pomeron with a relatively large intercept Ajp = ajp — 1 = 0.2 — 0.4, and exceedingly small slope 
a'jp ~ 0.02 GeV~~ 2 . A Pomeron with these characteristics appears in N=4 SYM [6-10] with a large coupling, 
which, at the moment, is the only theory that allows us to treat the strong interaction on a theoretical basis. 
Having a'jp — > 0, the Pomeron in these models has a natural matching with the hard Pomeron that evolves 
from perturbative QCD. Hence, these models could be a first step in building a self-consistent theoretical 
description of soft interactions at high energy, in spite, of the many phenomenological parameters (of the 
order of 10-15) appearing in these models. 

The new LHC data on soft interaction scattering at high energy (see Refs. [11-14]) shows that this 
hope was premature. Indeed, all of the models [1-5] failed to describe the LHC data. At first sight, this 
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suggests the need for a comprehensive revision of the main ingredients of the Pomeron models applied to 
high energy soft interactions. 

However, an alternate approach is to adjust the Pomeron model's parameters. Our first attempt 
to do this was published in Ref. [15], In which we concluded: "In spite of the fact that the values of 
the parameters, extracted from our current fitting, are slightly different from our previous values, the 
overall picture remains unchanged. Our updated total and elastic cross sections are slightly lower than the 
published TOTEM values [14], but still within the relatively large experimental error bars. Should future 
LHC measurements confirm the present TOTEM values, we will need to revise our dynamic picture for 
soft scattering." 

In the present paper we retract the above statement, as we find that the above conclusion was prema- 
ture. The set of parameters in our previous paper was found by fitting all data with energy W > 500GeV, 
including the LHC data. In the present version of our GLM model we made no changes other than tuning 
the 3 Pomeron parameters. Our tuned Ajp changes from 0.21 to 0.23, while G^ip and 7, the Pomeron- 
proton vertices, are unchanged. The small change in the value of Ajp is sufficient to produce the desired 
results in our a to t and o e \ output values for LHC energies, while the changes in the output values of the 
other observables are small enough not to spoil the good reproduction of the data achieved in Ref. [15]. 

Relying on our model that describes all available data, we also include in this paper predictions for 
LHC higher energies (W = 8TeV and W = 14TeV). As there has been a recent renewal of interest in 
cross sections eminating from cosmic ray data [22], we also list our results for W = 57TeV. 

2. Our Model 

To make our presentation self-contained we give a brief overview of our model, in spite of the fact, that all 
the ingredients and formulae have been published (see Ref. [1,2, 15]). 

2.1 Two channel model 

To account for diffraction dissociation in the states with masses that are much smaller than the initial 
energy, we use the simple two channel Good- Walker model. In this model we introduce two eigenwave 
functions, ipi and ip2, which diagonalize the 2x2 interaction matrix T, 



The two observed states are an hadron whose wave function we denote by iph, and a diffractive state 
with a wave function i/jd. These two observed states can be written in the form 



A,k =< A ^fc|T|Vv ipk' >= A,k $i,i> Sk,k'- 



(2.1) 



ip h = a ipi + (3 i>2 , ipD = ~/3 ipi + a V>2 



(2.2) 



where, a 2 + ft 



'■2 



= 1 
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2.2 Eikonal approach 

Using Eq. (|2.1|), we can rewrite the s-channel unitarity constraints in the form 

2 Im A i>k (s, b) = \A i>k (s, b)\ 2 + G™ k {s, b), (2.3) 

where, Gfi denotes the contribution of all non GW inelastic processes. 
In a general solution of Eq. (2.3) 

A hk (s,b) =i(l-exp (_^lR)) , ( 2 .4) 



2 

Gf fe ( S ,6) = l-exp(-O ijfe ( S ,fe)), (2.5) 

in which fi,- ^ are arbitrary. In the eikonal approximation Qi jk are taken as being real. In our model we 
choose fij.fc as the contribution of a single Pomeron exchange (see Eq. (|2.7|)). 



From Eq. (2.5) we deduce that the probability that the initial state (i, k) remains intact during the 
interaction, is Pf k = exp (— Q itk (s, b)). 

It should be stress that in this model we describe the variety of states produced with different masses 
and other quantum numbers, in diffraction dissociation, by a single state (Vt>) with unknown mass. Using 



Eq. (|2.l| ),Eq. ( |2.2[) and Eq. (2/7) we obtain that one Pomeron exchage has the following contributions to 
the different processes: 



Elastic scattering: p + p — )> p + p A PiP = a A Ai t i + 2a 2 (3 2 Ai >2 + /? 4 ^-2,2; 

Single diffraction: p + p^ D + p A P)D = aj3 (-a 2 A lil + (a 2 - (3 2 )A lj2 + /3 2 A 2 ,2) \ 

Double diffraction: p + p -> D + D Ae>,d = ol 2 /3 2 (,4i,i - 2Ai, 2 + A 22 )- (2.6) 



Eq. (2.6) shows the particular realization of the Good- Walker mechanism [16] for diffraction production, 
in the two channel model. It is worthwhile mentioning that this mechanism is the only source of diffraction 
in N=4 SYM, due to AdS-CFT correspondence (see for example Ref. [10]). 

2.3 Pomeron 

As we have mentioned, our basic ingredient is the soft Pomeron, which has the following contribution to 
the elastic scattering amplitude for the process: i + k — >• i + k where i and k are the states that are described 
by the wave functions i^i and ipk 

tt i)k (s,b) = gi (b)g k (b)P(s). (2.7) 
P(s) = s A , and gi(b) and g k (b) are the Pomeron-hadron vertices parameterized in the form: 

gi (b) = 9l St (b) = ^-mfbKx (m,6) . (2.8) 
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Si(b) is the Fourier transform of ( 1+g 2/ m i)2 ; where, q is the transverse momentum carried by the Pomeron, 
I. I.K. 

The form of P(s) that we use, corresponds to a Pomeron trajectory slope a'jp = 0. This is compatible 
with the exceedingly small fitted value of a'jp, and in accord with N=4 SYM [6]. 

fij^ is the imaginary part of the scattering amplitude for a single Pomeron exchange. 
2.4 Pomeron interactions 

In the case of Ajp — > 0, we know that the Pomeron interaction leads to a new source of diffraction 
production with large mass (M oc s), which cannot be described by the Good- Walker mechanism. For 
Ajp > 0, the Pomeron interaction is responsible for diffraction production of finite (independent of the total 
energy) mass: log (M 2 /so) °c 1/Ajp. This source can be treated in the same manner as the Good- Walker 
mechanism (see Ref. [17]). We believe that we need to account for the Pomeron interaction separately, 
since the typical diffractive mass depends on the Pomeron intercept, and this effect needs to be taken into 
account in the fitting procedure. 

Taking a'jp = 0, allows us to sum all diagrams having Pomeron interactions [1,2]. This is the great 
advantage of such an approach. In our model we only take into account triple Pomeron interaction vertices 
(G31P), this provides a natural matching to the hard Pomeron, since at short distances G^jp oc a|, while 
other vertices are much smaller. For a thorough description of the procedure for summing all diagrams, 
we refer to our papers (see Refs. [1,2,15]. 

Enhanced diagrams: In our model [2], the Pomeron's Green function that includes all enhanced dia- 
grams, is approximated using the MPSI procedure [18], in which a multi Pomeron interaction (taking into 
account only triple Pomeron vertices) is approximated by large Pomeron loops of rapidity size of Ins. 

We obtain 

Gjp (Y) = 1 - exp f,^) ^ r (0, -1—) , (2.9) 



T(Y) J T(Y) \ ' T(Y) / 
in which: 

T(Y) = 7 e AjpY . (2.10) 
r (0, l/T) is the incomplete gamma function (see formulae 8.35 in Ref. [19]). 

Semi-enhanced (net) diagrams: A brief glance at the values of the parameters of our model (see 
Ref. [15]), shows that we have a new small parameter, namely, G^jpP(s) <C 1, while G31P gi P(s) ~ 1. 

Gsip g% P(s) I , but do not contain any of G^p P(s) 
contributions, net diagrams. Summing the net diagrams [1], we obtain the following expression for fij fc(s, 6): 



rf*(Y;b) = d 2 b' V J - ' V . J , V J ' ~ - ^ "> . (2.11) 



9 i(b') g k (b-b>) (lhGjp(T(Y))) 
1 + (G3W7) Gip(t(Y)) \ 9i (v) +g k (b- V 

Gsjp is the triple Pomeron vertex, and 7 2 = / i^rG^jp. 
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2.5 Final formulae 

For the elastic amplitude we have: 

a el (b) = i (a 4 A 1A + 2a 2 p 2 A 1>2 + /3 4 A 2j2 ) . (2.12) 

For diffraction production we introduce an additional contribution due to the Pomeron interaction, which 
we call non-Good- Walker (non-GW). For single diffraction we have (see Ref. [15]): 

Ag kjL r = f d 2 b' 2 A -i) 9i (b - b', rm) 9l (b', m,) g k (b', m k ) 

x Q (gk,™>k$, Y -Y m ) Q (gi,mi,b',Y -Y m ) , (2.13) 

where, 

Q(g,m,b;Y) = - J^TT^IVS oTh V ( 2 - 14 ) 

1 + (G3JP/7) gGip{Y) S (b,m) 

For double diffraction we have: 

r~Ai% = f d 2 b'A (b-b',mi) g k (b',m k ) 

„2A SY . 



x Q [gi,mi, b - b,Y - Yi) e Q (g k ,m k , b',Yi - SY) . (2.15) 



Eq. ( 2.1 3| ) and Eq. ( 2.15| ) are the simplifications of the exact formulae of Ref. [1]. We have checked 



that they approach the values of the exact formulae reasonably well, within 5 — 10%. 

For single diffraction, Y = In (M 2 /so), where, M is the SD mass. For double diffraction, Y — Y\ = 
In (Mi /s ) and Yy-SY = In (M|/s ), where Mi and M 2 are the masses of the two bunches of hadrons 
produced in double diffraction, sq is the minimal produced mass, which is about 1 GeV . The integrated 
cross section of the SD channel is written as a sum of two terms: the GW term, which is equal to 

= j d 2 b(a(3{-aa 2 A l!l + (a 2 -/3 2 )A l! 2 + (3 2 A 2j2 }) 2 . (2.16) 

The second term describes diffraction production due to non-GW mechanism: 

afp W = 2 jdY m fd 2 b (2.17) 

{ a 6 Affi,! e -°M(^) + a 2 /? 4 < 2 ,2 e- n ^ Y ^ + 2 a 4 /3 2 A s ^ h2 g-^M^+^O^)) 
+ (3 2 a 4 A s 2 d 1:1 e- n ^ (Y ^ + 2 /?V A|f 1>2 e -^MY;b)+n 2 , 2 (Y;b)) + p e A ^ e -n 2 , 2 (Y ; 6) j 

The cross section of the double diffractive production is also a sum of the GW contribution, 

= / rf 2 &a 2 ^ 2 {^i,i - 2Ai i2 + ^ 2 , 2 } 2 , (2.18) 
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to which we add the term which is determined by the non GW contribution, 

4 GW = J d% [a 4 ^ ie -^, l( y;6) + 2a 2 p 2 A{% e - Q ^ Y ;b) + ^ A M e -n 2 , 2 (Y;b) j (219) 

In our model, the GW sector contributes to both low and high diffracted mass, while the non-GW sector 
contributes only to high mass diffraction (log (M 2 /so) ~ 1/Ajp). 

3. Results 

The output of our current model is presented in Fig. 1, in which we display our calculated cross sections 
and forward elastic slope. It is interesting to compare the quality of our present results with our previous 
output [15]. Recall that the two fits have almost identical values for the free parameters, with a single 
change of Ajp from 0.21 to 0.23. 
We list the main features of our results: 

• The main feature of our present calculation is the excellent reproduction of TOTEM's values for 
a tot and a e {. The quality of our good fit to B e i is maintained. As regards Oi„ e ;, our results are in 
accord with the higher values obtained by ALICE [11] and TOTEM [14]; ATLAS [12] and CMS [13] 
quote lower values with large extrapolation errors. We refer the reader to [25] who suggests that the 
lower values found by ATLAS and CMS maybe due to the simplified Monte Carlo that they used to 
estimate their diffractive background. 

• The quality of our output at lower energies, when compared with ISR data, is not as good as our 
previous results [3], but still acceptable. Recall that Reggeon exchange, which is included in our 
model, plays an important role at the low energy end of our data, and a negligible role at higher 
energies. As our main goal is to provide a good description of the LHC data, we have not tuned the 
Reggeon parameters, which could lead to an improved characterization of the ISR measurements. 

• An interesting observation is that our updated output strongly supports the CDF total and elastic 
cross sections rather than the E710 values. This is a common feature of other models [4,21], that 
have succeded in reproducing the TOTEM results by making a radical change in their modellings. 

• Note that our model is the only one which offers a good reproduction of ISR diffraction, and a 
reproduction of the diffractive cross sections at higher energies. 

• Our reproduction of SD and DD cross sections is complicated by the lack of common definitions of 
signatures and mass bounds on the diffractive components. All models on the market have intro- 
duced at least two different mechanisms to describe diffraction production. In our model these two 
mechanisms are : the Good- Walker production of the diffraction state with finite unspecified mass 
(independent of energy and of the values of the parameters in our model); and the diffraction due to 
Pomeron interactions where the typical mass depends on Ajp, and which we call non Good- Walker 
mechanism. In other models (see Refs. [3-5]) the two different mechanisms are called: low mass 
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diffraction and high mass diffraction. It is assumed that the diffractive production of mass less than 
Mo ~ 2 — 3 GeV stems from Good- Walker mechanism, while the production of masses larger than 
Mq is due to the Pomeron interaction. This partition appears natural in the case where Ajp — > 0, 
when the Pomeron interaction is responsible for diffractive production in the region where the mass 
increases with energy. However, for Ajp ^ both mechanisms lead to the production of final states 
with typical masses, and therefore the artificial requirement that the Good- Walker mechanism pro- 
duces masses less than Mq, appears unnatural. In our approach we do not specify the value of the 
mass in the Good- Walker mechanism, but it is certainly larger than Mq = 2 — 3 GeV. As the Pomeron 
interaction gives a negligible contribution at W = 900 GeV, we estimate that our typical mass for 
the GW mechanism is about 10 — 20 GeV. Such a large mass for the GW contribution, explains why 
our Pomeron interaction leads to a small contribution to the cross section for diffraction production 
due to Pomeron interactions. Note that the experimental errors on the diffractive cross sections are 
relatively large, and thus the data does not provide a good enough resolution to differentiate between 
competing models. 

• Our calculated da e \(t < 0.55GeV 2 )/dt is presented in Fig.2 together with the corresponding data. 
The quality of the fit is very good in this region of t. It shows that the impact parameter dependence 
of the model does not change considerably and reproduces the experimental data as well, as in our 
previous version of the model. We refrain from trying to reproduce the diffractive dip and higher t 
cross sections since our model is confined to the forward cone. 

• Table 1 summarizes our calculated cross sections and B e \ at 1.8 - 57TeV. The table provides the 
predicted values of the cross sections and forward slope at 8, 14 and 57 TeV. Note that the table 
indicates a slow approach of our predicted elastic amplitude to saturation. A common signature for 
saturation is a e \ = (Tinei = O.batot- Considering the energy behaviour of cr ine i/atot, the values of 
this ratio given by our model suggest a very slow increase to saturation, with (Tinel/ a tot = 0.77 at 
Tevatron energies decreasing to 0.73 at 57 TeV. 

The origin of such a slow approach to the black disc limit, turns out to be the same as with our previous 
set of parameters (see Ref. [15]): where one of the partial scattering amplitude, A\ t i (s,£>) < 1 at all 
energies higher than 1.8 TeV, while (s, b) ~ Aip. ( s > b) = 1 for b = at the same energies. 

• We note that our new set of parameters with Ajp = 0.23 does not change any qualitative features 
of the scattering amplitude. In particular, the impact parameter dependence in our model is similar 
to the b dependence of our previous model [15]. Therefore, it is not surprising that we obtain good 
agreement with the experimental data for da e i/dt as a function of t, at different energies (see Fig. ^). 

• The increase of Ajp leads to stronger shadowing corrections, that manifest themselves in smaller 
values for the survival probabilities for the large rapidity gap (S 2 ), compared to what we previously 
obtained with Ajp =0.21 [15]. The S 2 calculated with our present parameters are: 

S 2 = 7.02% [9.76%] atW = 1.8 TeV; S 2 = 2.98% [6.3%] atW = 7 TeV; 
S 2 = 2.7% [5.9%] at W = 8TeV; S 2 = 1.75% [3.66%] atW = 14 TeV; (3.1) 
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[. . . ] denotes the values obtained using our previous set of parameters (see Ref. [15]). 



W = V"* TeV 


1.8 


7 


8 


14 


57 


(Jtot mb 


79.2 


98.6 


101 


109 


130 


a e i mb 


18.5 


24.6 


25.2 


27.9 


34.8 


<r sd (M 2 < 0.05s)to6 


8.2 + [2.07] 


10.7 + [4.18] 


10.9 + [4.3] 


11.5 + [5.81] 


13 + [8.68] 


Odd mb 


5.12 + [0.38] 


6.2 + [1.24] 


6.32 + [1.29] 


6.78 + [1.59] 


7.95 + [5.19] 


B el GeV- 2 


17.4 


20.2 


20.4 


21.6 


24.6 


<7 ine l mb 


60.7 


74 


75.6 


81.1 


95.2 


%\ t=0 mb/GeV 2 


326.34 


506.4 


530.7 


608.11 


879.2 



Table 1: Predictions of our model for different energies W. The quantities appearing in square brackets [. . . ] denote 
the contribution due to Pomeron interactions (so called non-GW part of the cross section) 

4. Conclusions 

Our goal in this paper is to deliver the message: we have succeded in building a model for soft interactions 
at high energy, which provides a good description all high energy data, including the LHC measurements. 
This model is based on the Pomeron with a large intercept (Ajp = 0.23) and slope a'p « 0. We find no 
need to introduce two Pomerons: i.e. a soft and a hard one. The Pomeron in our model provides a natural 
matching with the hard Pomeron in processes that occur at short distances. 

The attractive feature of our model is that we are able to describe data over the whole energy range. 
Our description of the data at lower energies is slightly worse than our description of the data at high 
energy W > 900 GeV . In addition we note that the qualitative features of the model are close to what we 
expect from N=4 SYM, which is the only theory that is able to treat the long distance physics on a solid 
theoretical basis. 
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Figure 1: Comparison with the experimental data the energy behaviour of the total (Fig. |l|-a), elastic (Fig. |T|-b), 
single diffraction (Fig. |l|-c), double diffraction (Fig. Q-d) and inelastic (Fig. Q-f) cross sections and elastic slope( 
Fig. 0-e), The solid lines show our present fit . The data has been taken from Ref. [20] for energies less than the LHC 
energy. At the LHC energy for total and elastic cross section we use TOTEM data [14] and for single and double 
diffraction cross sections are taken from Ref. [If]. 
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Figure 2: da e i/dt versus \t\ at Tevatron (blue curve and data)) and LHC ( black curve and data) energies (W = 
1.8 TeV , 8TeV , 7Tev and 14 TeV respectively). Data from Refs. [14,23]. 
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